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COMPASS is a fixed-target experiment at CERN SPS which investigates the structure and spec- 
troscopy of hadrons. During nine weeks in 2008 and 2009, a 190GeV/c proton beam impinging 
on a liquid hydrogen target was used in order to study the production of exotic mesons and glue- 
ball candidates at central rapidities. As no bias on the production mechanism was introduced 
by the trigger system, the contribution from diffractive dissociation of the beam proton poses a 



challenge. We select a centrally produced sample by kinematic cuts and introduce a model to de- 



scribe the data in terms of partial waves. Preliminary fits are presented, which are consistent with 
results from previous experiments. Particular attention is paid to the ambiguities in the amplitude 
analysis of the two-pseudoscalar final state. 
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1. Introduction 

Quantum Chromodynamics predicts objects composed entirely of valence gluons, so-called 
glueballs. Their existence, however, could not be confirmed experimentally to the present day. 
One of the goals of the COMPASS experiment [1] at CERN is to study the existence and signatures 
of glueballs, in continuation of the efforts that were made at the CERN Omega spectrometer in the 
late 1990s [2]. Since Pomerons are considered to have no valence quark contribution, Pomeron- 
Pomeron fusion was proposed to be well suited for the production of glueballs. This process can 
be realised in a fixed-target experiment by the scattering of a proton beam on a proton target, where 
a system of particles is produced centrally (cf. Figure 1). 



Approximately 30% of the total COMPASS data recorded with a 190GeV/c proton beam 
impinging on a liquid hydrogen target was used for the presented analysis. In order to separate the 
central 7r + 7r~ system from the fast proton pf, a cut on the invariant mass combinations M(pn) > 
1.5GeV/c 2 was introduced. The effect of this kinematic cut on the Feynman xf distributions for 
the fast (yellow) and the slow proton (red) as well as the di-pion system (blue) is shown in Figure 2. 
The % + %~ system lies within \xf\ < 0.25 and can therefore be considered as centrally produced. 

However, central production comprises several production processes, their characteristic sig- 
nature being the different dependence of the cross section on the centre-of-mass energy y/s of the 
reaction. While Pomeron-Pomeron scattering should be only weakly dependent on s, Pomeron- 
Reggeon scattering is predicted to scale with 1 / y/s and Reggeon-Reggeon scattering with 1 /s [3]. 
This behaviour can be observed experimentally, e.g. by comparing the mass spectra obtained at 
COMPASS energy to those the Omega spectrometer [4] obtained at different energies (cf. Fig- 
ure 3). 

The dominant features of the n + K~ invariant mass distribution, the p(770), the /2(1270) and 
the sharp drop in intensity in the vicinity of the /o(980), can be observed at all three different 
centre-of-mass energies. The relative yield of p-production decreased rapidly with increasing y/s, 
which is explained by a Reggeised pion-pion contribution (« l/s 2 ). On the the other hand, the 
enhancement at low masses as well as the /b(980) remain practically unchanged, a fact that is 
characteristic for ^-independent Pomeron-Pomeron scattering. 
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Figure 1: Central production 



Figure 2: Feynman Xf distributions 
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Figure 3: Invariant mass of the % + % system for three beam energies 



2. Partial-Wave Analysis 

The partial-wave analysis has been performed assuming that the % + %~ system is produced 
by the collision of two particles emitted by the scattered protons. They are also referred to as 
exchanged particles, and carry the squared four-momentum transfer t\ from the beam proton and tj 
from the target proton, respectively. 

We make the strong assumption that t\ only transmits the helicity X = 0. In this limit, t\ can 
be seen as a vacuum-like Pomeron with J PC = ++ , which is treated like an external particle. The 
beam and the fast outgoing proton pf are irrelevant in this picture. If we accept the space-like 
Pomeron as the incoming beam, we can construct a Gottfried- Jackson frame [5] for the di-pion 
system. Hence, the exchanged particle t\ in the centre-of-mass frame of the 7t + 7t~ system defines 
the z-axis for the reaction. We have to emphasise that we fixed the choice to t\ in contrast to 
previous experiments (e.g. [6]) to be able to correct for the different t acceptances of the trigger. 

The v-axis of the right-handed coordinate system is defined by the cross product of the mo- 
mentum vectors of the two exchange particles in the pp centre-of-mass. Two variables specify the 
decay process, namely the polar and azimuthal angles (cos(0), 0) of the ii~ in the di-pion centre- 
of-mass frame relative to these axes. Figure 4 illustrates the distribution of these decay variables as 
a function of the di-pion mass. 

The decay amplitudes for a two-pseudoscalar system with relative angular momentum i and 
its projection on the quantisation axis m is given by the spherical harmonics Yf n (6, (j)) in spherical 
coordinates. To profit from the fact that the strong interaction conserves parity, a reflection operator 
was introduced in [7] as a parity operator followed by a rotation around the j-axis such that all 
momenta relevant to the production process are recovered. Introducing a quantum number e which 
is called reflectivity and which can have the values ± 1 , the eigenstates of this reflection operator 
can be constructed as: 



with a normalisation constant c m . The two classes of states s = ±1 correspond to different produc- 
tion processes in the asymptotic limit [5] and can consequently not interfere. 



Y^(6,<l>) = c m Y l m 



(6,<l>)-s(-l) m Y[ m (6,<l>) 



(2.1) 



3 



Partial- Wave Analysis of the Centrally Produced 7T + 7T System in pp Reactions at COMPASS 
Tobias Schliiter 



o- 

-0.2 ~— 
-0.4 r 
-0.6 ~— 
-0.8 — 



3.2 0.4 0.6 0.8 



.2 1.4 1.6 1.8 2 
Invariant Mass of %*%' (GeV/c 2 ) 




0.2 0.4 0.6 0.8 



Invariant Mass of Jt + Jt" (GeVAr) 



(a) cos(0) (b) 

Figure 4: Decay variables as a function of the K + n~ mass 



With the complex transition amplitudes T e i m and an explicit incoherent sum over the reflectiv- 
ities, the intensity in mass bins can be expanded in terms of partial waves: 



/(M) = E 



(2.2) 



We adopt the notation 7^ from [6] to construct the basic wave set ,Pq,P± ,Dq ,D~[ and 
P^jD^. Since the overall phase for each reflectivity is indeterminate, one wave in each class can 
be defined real, which means the imaginary part of that transition amplitude is fixed to zero. These 
so-called anchor waves are emphasised in bold font. 

An extended likelihood fit in 20MeV/c 2 mass bins is used to find the parameters T e i, n , such 
that the acceptance corrected model I(6,<p) matches the measured data best. 



3. Ambiguities 

As shown in [8, 9], the angular dependency of the intensity I(Q.) for two-pseudoscalar final 
states can generally be expressed in terms of \G{u) | 2 by the introduction of a variable u = tan(0 /2). 
The function G(u) is a complex polynomial of the order of 11, where I is the the highest considered 
spin in the system. This polynomial can be factorised in terms of its complex roots u^, the so-called 
Barrelet-zeros [10]. Since the function G(u) only enters as absolute square in the expression for 
the angular distribution, the complex conjugate of a root u* k is an equally valid solution. That 
means, that there are in general 2 2e ~ l different ambiguous solutions which result in exactly the 
same angular distribution. This ambiguity has to be resolved by physical arguments. 

The system of S, P and D waves used in the presented analysis has eight ambiguous solutions. 
By using different random starting values for the fit to the same data and MC sample, it was shown 
experimentally that very different solutions can be obtained. However, the fitted production am- 
plitudes r e / m for one single attempt can be used to calculate all eight solutions analytically [9] . In 
order to find the complex polynomial roots numerically, Laguerre 's method was applied. Figure 5 
illustrates the real and imaginary parts of these four roots for all mass bins, where a sorting depend- 
ing on the real part has been performed. They are well separated from each other and do not require 
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Figure 5: The (a) real and (b) imaginary parts of the Barrelet-zeros as a function of the % + % mass 



a special linking procedure. The imaginary parts do not cross the real axis, hence bifurcation of the 
solutions does not pose a problem and the solutions can be uniquely identified. 

By fixing u\ and allowing 1*2,3,4 to undergo complex conjugation, the entire set of eight am- 
biguous solutions is calculated. For six solutions, most of the intensity is formed by one single 
wave. This is clearly unphysical, since we know of at least one resonance in each of the spins 0, 1 
and 2. The choice among the remaining solutions is not evident, the intensity and phase distribu- 
tions are very similar. 
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Figure 6: Solution compatible with physical constraints, intensities (blue) and relative phases (red) 
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As a last step, the calculated values for one chosen solution are reintroduced to the fit as starting 
values, which probes the convergence and provides the correct covariance matrix. Figure 6 shows 
the intensities and phases for one solution compatible with the physical constraints, depicted in a 
matrix form. A clear signal from /2(1270) can be seen in the Dq wave, accompanied by a phase 
motion with respect to the 5q . There is also evidence for p(770) in both and Pf. At masses well 
below 1 GeV/c 2 , however, the scalar wave does not show the expected threshold enhancement. 

To evaluate the fit quality, the decay amplitudes of phase-space Monte-Carlo events are weighted 
by the production amplitudes 7} obtained in the minimisation. The projected angular distributions 
show very good agreement with the data. 

4. Outlook 

With these proceedings, we want to show that we are able to select a centrally produced sample 
and describe the main features of the data in terms of partial waves. Using the methods from [6], 
we can reproduce the analysis with comparable results. 

The complicated problem of central production was approached following a model employed 
by WA102, an earlier experiment at CERN, inheriting some limitations and assumptions. In the 
meantime, other ideas were developed based on their data and analyses which can help to determine 
spin and parity of particles produced in central exclusive processes [11]. In addition to measuring 
the decay products, it can have particular advantages to study angular correlations between the 
outgoing protons. In the long run, we might study the data with a combination of both. 

In order to obtain a complete picture of the light scalar meson sector, the decay of identified 
resonances into different final states is needed to understand their composition. Since the trigger 
was not selecting a particular final state and the generalisation of the analysis of the two-pion 
system to other two-pseudoscalar final states is trivial, the next natural step will be the analysis of 
the K + K~ system. Other available channels are K S K S , 7i Q K° and t\t\. We expect that COMPASS 
will be able to shed light on this sector, outperforming previous experiments in terms of data as 
well as resolution. 
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